Spin qubits in double quantum dots - entanglement versus the Kondo effect 
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We investigate the competition between pair entanglement of two spin qubits in double quantum dots attached 
to leads with various topologies and the separate entanglement of each spin with nearby electrodes. Universal 
behavior of entanglement is demonstrated in dependence on the mutual interactions between the spin qubits, the 
coupling to their environment, temperature and magnetic field. As a consequence of quantum phase transition an 
abrupt switch between fully entangled and unentangled states takes place when the dots are coupled in parallel. 
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Introduction.- After the recent discovery of quantum com- 
puting algorithms, their practical potential led to the inter- 
est in quantum entanglement spurred on by the fact that if a 
quantum computer were built, it would be capable of tasks 
impracticable in classical computing 1 1]. Nanostructures con- 
sisting of coupled quantum dots are candidates for required 
scalable solid state arrays of electron spin qubitsQ, 0. The 
interaction of such qubits with the environment is in general a 
complicated many-body process and its understanding is cru- 
cial for experimental solid state realisation of qubits in single 
and double quantum dots (DQD)|4]. Recent experiments on 
semiconductor double quantum dot devices have shown that 
electron occupation may be controlled down to the single- 
electron level by surface gates|5]. Also spin entangled states 
were detected|6], DQDs were used to implement two-electron 
spin entanglement|7], and coherent manipulation and projec- 
tive readout|8] was demonstrated. 

The purpose of entangled qubit pairs is to convey quantum 
information through a computing device|l]. The entangle- 
ment of two spin qubits may be uniquely defined through von 
Neuman entropy or, equivalently, concurrence patr 
of qubits may be realised, e.g., as two separate regions, each 
occupied by one electron in a state |s)a,b of either spin, s =| 
or j. For a system in a pure state |*ab) = J2 SS ' a ss' I s ) A ® 
|s')b, the concurrence as a quantitative measure for (spin) en- 
tanglement is given by| 10] Co — 2\a^a^ — a^an\. Two 
qubits are completely entangled, Co = 1, if they are in one of 
the Bell states®, e.g., singlet |* A b) oc | fJ,) - | IT}- 

The setup and main results- We focus on entanglement be- 
tween two electrons confined in two adjacent quantum dots 
weakly coupled by electron tunneling in a controllable man- 
ner, Fig.^a). The inter-dot tunneling matrix element t deter- 
mines not only the tunneling rate, but also the effective mag- 
netic superexchange interaction J ~ 4t 2 /U, where U is the 
scale of Coulomb interaction between two electrons confined 
on the same dot. By adjusting a global back-gate voltage, ex- 
actly two electrons can be confined to the dots A-B on aver- 
age. 

Additional gate voltages are applied to independently con- 
trol tunneling to the electrodes, t n . Depending on the val- 
ues of t n , various topologies can be realised and even for 
very weak coupling, the spin of confined electrons may be 




Figure 1 : (color online) (a) Double quantum dot system, t and t n are 
the matrix elements for tunneling between dots A and B and from 
the dots to the electrodes, respectively, (b) Entanglement between 
two spins on the quantum dots as a function of the interdot exchange 
coupling J: below J c , the two spins are not entangled and the DQD 
is in some type of the Kondo regime. For elevated temperatures and 
magnetic field above J c , the entanglement is zero if J < max(T, B). 



screened due to the Kondo effect, where at temperatures be- 
low the Kondo temperature Tjc a spin-singlet state is formed 
between a confined electron and conduction electrons close to 
the Fermi energy. Conductance and some other properties of 
such systems have already been studied, without considering, 
however, the analysis of the entanglement and its relationship 
to the many-body phenomena embodied in the Kondo effect. 

Qualitatively the physics related to qubit pairs in coupled 
DQDs is plausible and can be summarised as: 

(i) If t/U is not small the electrons tunnel between the dots 
and charge fluctuations introduce additional states with zero 
or double occupancy of individual dots JTTl IT2I1 . Due to sig- 
nificant local charge fluctuations this regime is not particularly 
appropriate for the spin-qubits manipulation. 

(ii) For systems with strong electron-electron repulsion, 
charge fluctuations are suppressed and the states with single 
occupancy - the spin-qubits - dominate. Due to the effective 
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antiferromagnetic Heisenberg interaction the spins A and B 
tend to form the singlet state. The physics of such qubit pairs 
may be compared to the two-impurity Kondo problem studied 
by Jones, Varma and Wilkins two decades ago| 13]. There two 
impurities form either two Kondo singlets with delocalised 
electrons or bind into a local spin-singlet state which is vir- 
tually decoupled from delocalised electrons. The crossover 
between the two regimes is determined by the relative values 
of the exchange energy J and twice the Kondo condensation 
energy, of order the Kondo temperature Tk- 

As shown in this paper, our numerical results for represen- 
tative DQD systems with different topology of possible ex- 
perimental realisations reveal much more diverse physical be- 
haviour. However, (i) in all cases the spin qubits are unen- 
tangled for J below some critical value J c , where the actual 
value of J c crucially depends on the setup topology, and (ii) 
at elevated temperatures T > and external magnetic field 
B ^ the entanglement is additionally suppressed and gener- 
ically zero when J < max( J c , T, B), as schematically shown 
in Fig. Eft). 

Quantitative results.- For simplicity we model DQD us- 
ing the two-site Hubbard Hamiltonian H = —t X) s ( c As c Bs + 
Cq s cas) + U Xyj=A b n i\ n ili where c\ s creates an electron 
with spin s in the dot i = A or i = B and n, s = c\ s Ci S is 
the number operator. The dots are coupled to the left and right 
noninteracting lead as shown in Fig.^a). 

DQDs as considered here can not be described with a 
pure quantum state and concurrence is not directly given by 
Co- It is related to the reduced density matrix of the DQD 
subsysteml 10Hl4lll5ll . where for systems that are axially sym- 
metric in spin space the concurrence may conveniently be 
given in the closed formlUfill. 

C = max(0,C T ;,C||)/(P n +P||), (1) 

C n = 2\(S+S B )\-2^f(P{pl)(Pipt}, 

q, = 2|(5+5+)|-2 V /(pTpi}(p|p^, 

where = (S~y = c|j c a is the electron spin raising op- 
erator for dot i = A or B and P? = nj s (l — n^-g) is the 
projection operator onto the subspace where dot i is occupied 
by one electron with spin s. P^ = (P a Pb + P\Pb) and 
P|l = (P A Pj3> + P\Pb) are probabilities for antiparallel and 
parallel spin alignment, respectively. 

We have determined concurrence for all three possible 
topologically non-equivalent two terminal experimental ar- 
rangements: double quantum dots (i) coupled in series, (ii) 
laterally side coupled, and (iii) coupled in parallel. Concur- 
rence was determined numerically from Eq. Q, where expec- 
tation values correspond to a many-body state with the chem- 
ical potential in the middle of the electron band, which guar- 
antees that the dots are singly occupied, (?ia,b) = 1. Our ex- 
tensive investigation over the full parameter range for various 
topologies indicates that all show generic behaviour outlined 
in Fig. but quantitatively can differ by many orders of 



magnitude, which should be taken into consideration in exper- 
iments with such DQD. Numerical methods were based on the 
Gunnarsson-Schonhammer (GS) proiection-operatorl 17 1 1 1 3. 
ll% and numerical renormalisation group (NRG) 12(1 l2ll 12211 
methods. 

Serially coupled DQD.- First we consider serially coupled 
DQD, which models entangled pairs that may be extracted us- 
ing a single-electron turnstile|23]. Here — t' and <2,3 = 
with the hybridisation width of each dot T = (i') 2 /to, where 
4io is the bandwidth of noninteracting leads. Entanglement of 
a qubit pair represented with quantum dots in the contact with 
the leads (fermionic bath) was not quantitatively determined 
so far, although this system has already been extensively stud- 
ied 1 18, 24, 25] (and references therein). 

In analogy with entanglement at zero temperature studied 
recently in a many-body ground state|2g], we consider here 
concurrence of DQD at fixed temperature and static magnetic 
field along the z-axis. Expectation values (...) in the concur- 
rence formula Eq. Q correspond to thermal equilibrium of 
the system, therefore (S^Sq) — here. Qualitatively, the 
concurrence is significant when enhanced spin-spin correla- 
tions indicate inter-dot singlet formation. As shown in Fig. [2] 
the correlator (Sa ■ Sb) tends to —3/4 for J large enough to 
suppress the formation of Kondo singlets, but still J/U <C 1, 
that local charge fluctuations An A are sufficiently suppressed 
and Pfi + P|j — > 1. Concurrence, calculated for various val- 
ues of the Coulomb interaction strengths and in the absence 
of magnetic field is presented in Fig.|2j left bottom panel. As 
discussed above, C is zero for J < J\ c due to the Kondo ef- 
fect, which leads to entanglement between localised and con- 
ducting electrons 1 27] instead of the A-B qubit pair entangle- 
ment. In finite magnetic field irrespectively of temperature 
the concurrence abruptly tends to zero for B > J (not shown 
here)!!!. 

In particular, the local dot-dot singlet is formed and C > 
whenever singlet-triplet splitting J > J% c ~ 2.5Tk(T), 
where the Kondo temperature is given by the Haldane for- 
mula T K (T) = v / !7r/2exp(-7rC//8r). This is presented in 
the phase diagram in the (U /V, J/Tk) plane, Fig. [3] Dashed 
region corresponds to the regime of zero concurrence and is 
delimited by the line of the critical J = J\ c (red line). The 
charge fluctuations (Fig. [5] contour plot) are suppressed for 
sufficiently large repulsion, i.e., U/T > 10. In this limit 
and in vanishing magnetic field, the DQD can be described 
in terms of the Werner states|29] and becomes similar to re- 
cently studied problem of entanglement of two Kondo spin 
impurities embedded in a conduction band l!30ll . In this case, 
C n ~ 2(-(S A • S B ) - \) ~ P U ~ 2P|| for C n > 0. 
For large U/T, where the charge fluctuations vanish, the 
(Sa • Sb) = — | boundary (Fig. [3] dotted line) progressively 
merges with the C = line. 

Side-coupled DQD - In the side-coupled DQD configura- 
tion, t% t 2 = t', t3 t 4 = (Fig. 13 middle), the dot A is in 
direct contact with the electrodes, while the dot B couples to 
the conduction band only indirectly through the dot A. Be- 
cause the two electrodes are in contact only with dot A, is 
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Figure 2: (color online) Top panels: Spin-spin correlation function Si ■ S2, charge fluctuations on one quantum dot, An%, and probabilities 
Pf X> P\\ f° r P ^ Tk, P — > and with t! — to/ v20. Bottom left panel: concurrence C, corresponding to serially coupled dots, for a range of 
interactions U/T = 40, 32, 24, 16, 8, 4 and calculated with both, the NRG and the GS method (bullets), yielding the same Ji c , but due to the 
limited span of variational basis the GS method progressively overestimates C for U/T > 20 regime. Bottom right two panels: The results 
for side coupled and parallel configuration obtained from the NRG method. Temperatures range from the scale of the Coulomb repulsion 
parameter U, T/U — 0.4, to temperatures below the Kondo scale Tk', each consecutive curve corresponds to a temperature lowered by a 
factor 4. 
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Figure 3: (color online) In the phase diagram (U/T, J/Tk) full line 
separates C > and C = regions (line shaded), together with 
dotted line indicating (S A • S B ) = —1/4 (T — > and P = 0). 
Both lines merge for U/T > 12, where charge fluctuations An^ 
progressively become negligible (contour plot). 

r = 2(i') 2 /io, i-e- twice as much as in the previous case. 
Since Tk oc exp(— 7rC//8r), the Kondo temperature on dot A 
is strongly enhanced. 

For J > J 2c , the spins bind in an antiferromagnetic singlet, 
as in all other cases. For J < J 2c , the system enters the 'two 
stage Kondo' regime, characterised by consecutive screening 
of local moments 12 ll l3ll 13211 . At the Kondo temperature 
— Tk(T) the spin on the dot A is screened, while the 
spin on the dot B is compensated at a reduced temperature 

T [ k ] = exp(d 2 J/T^ ] ), (2) 



where are constants of order unity. The Kondo effect on 
dot A leads to the formation of a local Fermi liquid for temper- 
atures below The quasiparticle excitations of this Fermi 
liquid then participate in the Kondo effect on dot B at much 
lower Kondo temperature T K 1 3 1 ] . Such description is valid 
only when the temperature scales T^ 1 and tB' are widely 
separated. This no longer holds when J becomes compara- 
ble to T£\ see Eq. 0. The critical J 2c is thus still given by 
J 2c ~ Tjg. The crossover is very smooth and the transition 
from the inter-impurity singlet phase to the Kondo phase does 
not exhibit any sharp features. In fact, the low temperature 
fixed point is the same for J < J 2c and J > J 2c , unlike in the 
case of DQD in series. In the latter case, the Kondo phase and 
the inter-impurity singlet phase are qualitatively different and 
are characterised by different electron scattering phase shifts. 

In the J > J 2c singlet region, when the temperature is 
above J, the exchange interaction is too weak to bind the spins 
into a singlet and the entanglement is lost (see Fig. [2] bottom 
middle panel). In the J < J 2c Kondo region, the concur- 
rence is zero irrespective of temperature: for T < Tk it is 
zero due to the Kondo effect, and for T > Tk the spin-singlet 
cannot be restored, since T > J. Elevated temperature and 
magnetic field dependence is similar to the previous case of 
serially coupled dots 1 28]. 

Parallelly coupled DQD - In the case of parallel quantum 
dots (t n = t' and T = 2(t') 2 /to) the physics is markedly 
different from the case of the previous two configurations. 
The conduction band mediated effective Ruderman-Kittel- 
Kasuya-Yoshida (RKKY) interaction between the dots is, in 
our simplified model, ferromagnetic 12211 . Here the spins or- 
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der ferromagnetically into a triplet state (Fig. [2] right panels) 
and undergo a 5 — 1 Kondo screening at low temperatures in 
the regime 

64 r 2 

TX A U 

where c is a constant of order unity. This yields an uncom- 
pensated 5=1/2 residual spin and since half a unit of spin is 
quenched by the conduction band via the Kondo effect, there 
clearly cannot be any entanglement between the electrons on 
the dots. For J > |Jrkky|> the antiferromagnetic ordering 
wins and, once again, the electrons form an entangled sin- 
glet state. The transition from the Kondo phase to the singlet 
phase is in this case a true quantum phase transition !^ , and 
the concurrence drops abruptly from (surprisingly) C « 1, to 
C = when the exchange J between the dots is decreased 
below some Js c . Due to the sensitivity of concurrence on 
J ~ Jic can therefore parallely coupled DQD be considered a 
perfect entanglement switch. Critical coupling Js c is not de- 
termined by Tk, as previously, but rather by | JrkkyI- For 
this reason, concurrence drops to zero at a much higher tem- 
perature (compare middle and right bottom panels in Fig.|2ji, 
and is strictly zero for J < T as in previous two cases. In 
finite magnetic field C = if J < | Jrkky | + B (not shown 
here)!!!. 

If the couplings t n are not strictly equal, another Kondo 
screening stage may occur at low temperatures, in which the 
residual S — 1/2 spin is finally screened to zero|33]. In this 
case the quantum phase transition is replaced by a cross-over 
that becomes smoother as the degree of asymmetry between 
the couplings t n increases (results not shown). 

Conclussions.- We have found generic behaviour of spin- 
entanglement of an electron pair in double quantum dots. On 
the one hand, we have shown quantitatively that making the 
spin-spin exchange coupling J large by increasing tunneling 
i, leads to enhanced charge fluctuations, whilst on the other, 
a small interaction J < J c suppresses entanglement as the 
DQD system undergoes some form of the Kondo effect. Var- 
ious regimes are explained analytically and supported with 
typical numerical examples. In the limiting cases we found 

(i) two separate Kondo effects for serially coupled DQD; 

(ii) two-stage Kondo effect in side-coupled DQD; and (iii) 
5 = 1 Kondo effect with underscreening for parallelly cou- 
pled DQD, eventually followed by another 5=1/2 Kondo 
effect at lower temperatures. For two terminal setups, these 
are the only possible types of the Kondo effect; in a generic 
situation with arbitrary t n , one of these possibilities must oc- 
cur. 

In all cases, in spite of different Kondo mechanisms, the 
temperature and magnetic field dependence of entanglement 
is proven to be determined solely by the exchange scale J and 
not by the much lower scale of the Kondo temperature, which 
explains the universal behaviour of the entanglement shown in 
Fig-Eb). Critical J c , however, will for various experimental 
setups vary for several orders of magnitude. 
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